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Abstract 


The UVS instrument on the Juno mission recorded transient bright emission from 
a point source in Jupiter’s atmosphere. The spectrum shows that the emission is con- 
sistent with a 9600-K blackbody located 225 km above the 1-bar level and the duration 
of the emission was between 17 ms and 150 s. These characteristics are consistent with 
a bolide in Jupiter’s atmosphere. Based on the energy emitted, we estimate that the im- 
pactor had a mass of 250-5000 kg, which corresponds to a diameter of 1-4 m. By con- 
sidering all observations made with Juno UVS over the first 27 perijoves of the mission, 
we estimate an impact flux rate of 24,000 per year for impactors with masses greater than 


250-5000 kg. 


Plain Language Summary 


UVS is an ultraviolet spectrograph on NASA’s Juno spacecraft, which has been in 
orbit around Jupiter since 2016. In April 2020, UVS observed short-lived bright emis- 
sion from a point source in Jupiter’s atmosphere. The emission was consistent with a 
blackbody of temperature 9600 K. We suggest that this was a fireball produced by a 250- 


5000 kg meteoroid entering Jupiter’s atmosphere. 


1 Introduction 


As the largest and most massive planet in the Solar System, Jupiter undergoes a 
heavy bombardment of objects, ranging from tiny dust grains to kilometer-sized comets. 
The largest impacts, such as comet Shoemaker-Levy 9 (Asphaug & Benz, 1996) and the 
2009 impactor (Sdnchez-Lavega et al., 2010), occur rarely but leave scars on the planet 
that can persist for several months (Sdnchez-Lavega et al., 1998, 2011) and can affect 
the distribution of trace species in Jupiter’s upper atmosphere for decades (Lellouch et 
al., 2006; Cavalié et al., 2013; Benmahi et al., 2020). Impacts from objects in the 5-20 
m diameter range occur more frequently and produce short bright flashes of light that 
can be observed by amateur astronomers on Earth, but no observable debris (Hueso et 
al., 2018). Over a period of 8 years (2010-2017), amateur astronomers observed 5 of these 
smaller impacts, leading to an estimated impact rate of 10-65 per year (Hueso et al., 2018), 
an estimate that was not significantly modified by the observation of a 6th impact in 2019 
(Sankar et al., 2020). On the smallest scales, the constant influx of dust particles equates 


to hundreds of thousands of tons of material per year (Sreméevié et al., 2005; Poppe, 2016). 


While larger (>5 m) impacts can be observed from Earth, observations from or- 


biting spacecraft have the advantage of being able to detect the fainter flashes that are 


produced from the more frequent smaller impacts. One example of this is the small fire- 
ball observed by the camera on the Voyager 1 spacecraft, which Cook and Duxbury (1981) 
estimated was caused by an 11 kg (<0.5 m) meteoroid. In this paper, we present obser- 
vations of another fireball in Jupiter’s atmosphere, this time observed by the UVS in- 
strument on the Juno spacecraft. In Section 2, we describe the instrument and the method 
of observation. In Section 3 we discuss the properties of the bright blackbody-emission 
flash observed. In Section 4, we conclude that this flash was caused by a meteoroid of 
mass 250-5000 kg (diameter 1-4 m) entering Jupiter’s atmosphere, and we use the sum 


of all of our observations to estimate an impact flux rate. 


2 Observations 
2.1 Juno UVS 


The Ultraviolet Spectrograph (UVS) is an instrument on the Juno spacecraft, which 
has been in orbit around Jupiter since July 2016 (Bolton et al., 2017). UVS is a photon- 
counting imaging spectrograph, covering wavelengths of 68-210 nm in the far-ultraviolet 
(Gladstone, Persyn, et al., 2017). The main purpose of UVS is to study the morphol- 
ogy, brightness, and spectral characteristics of Jupiter’s auroras, and the instrument’s 
spectral range covers important H and Hg auroral emissions (Gladstone, Versteeg, et al., 


2017; Bonfond et al., 2017). 


Juno is a spin-stabilized spacecraft with a rotation period of 30 s; as the spacecraft 
rotates, the UVS instrument slit sweeps across Jupiter. The wavelength, slit position and 
precise time of UV photon detection events are recorded and this information is then used 
to build up spatial maps of the ultraviolet radiation (see Figure 1 for an example from 
a single spin). The UVS instrument slit consists of two wide segments on either side of 
a narrow segment. The wide parts of the slit have a width of 0.2° and a spectral reso- 
lution of 2.0-3.0 nm, while the narrow part of the slit has a width of ~0.025° and a spec- 
tral resolution of ~1.3 nm (Greathouse et al., 2013). The spatial resolution along the slit 
is 0.1-0.3°. The observations described in this paper were made at the top edge of the 


upper wide slit, where the spatial resolution is ~0.3°. 


2.2 Bright flash 


On occasion, the Juno UVS instrument has observed short-lived, localized ultra- 
violet emission outside of the auroral zone. In Giles et al. (2020), we described a set of 
eleven bright transient flashes that were observed with Juno UVS and shared similar char- 


acteristics; they lasted ~1.4 ms, they were dominated by Hə emission and the source re- 


Figure 1. A spatial map of Juno UVS measurements during a single spacecraft spin on 10 
April 2020. Photons with wavelengths 130-160 nm are shown in green and photons with wave- 
lengths 170-200 nm are shown in red. The inset image magnifies the area around the bright 


spot. 


gion was located ~260 km above the 1-bar level. Based on these characteristics, we con- 
cluded that these were likely to be elves, sprites or sprite halos, forms of Transient Lu- 
minous Events (TLEs) that occur in the upper atmosphere in response to tropospheric 


lightning. 


In this paper, we focus on another transient bright flash observed by UVS, but this 
observation has very different characteristics than the aforementioned TLEs. The ob- 
servation was made on 10 April 2020 at 12:57:10 UTC and the bright spot was located 
at a planetocentric latitude of 53°N and a System III longitude of 200° W. It was observed 
at an emission angle of 43° and a solar incidence angle of 125° (i.e. on the planet’s night- 


side). At the time of the observation, the spacecraft was at an altitude of 80,000 km above 


the 1-bar level of the planet and the sub-spacecraft latitude and longitude were 63°N and 
243°W respectively. This location was not observed by Juno’s infrared instrument (JIRAM, 
Adriani et al., 2008) or the visible light camera (JunoCam, Hansen et al., 2017) during 

this time period or afterwards. Figure 1 shows the UVS map from the spin in which the 
bright spot was recorded and includes a magnified image of the spot itself. This mag- 
nified image has dimensions 3° x3° on the sky, with each pixel being 0.1°x0.1°. The spot 
was only observed during a single spin; it was not seen two spins earlier or three spins 


later, which were the closest times when the same latitude and longitude was observed. 


The colors used in Figure 1 show the number of photon detections made in two dif- 
ferent spectral regions. Green is used to represent photons with wavelengths 130-160 nm 
and red is used to represent photons with wavelengths 170-200 nm. The UVS swath in- 
cludes a segment of Jupiter’s northern auroral oval and the auroral emission appears purely 
green in the image; this is because the main auroral Hə emission bands are at <170 nm. 

In contrast, the bright spot appears mostly yellow, indicating that there is significant 


emission at longer wavelengths as well. 


Although this bright spot is located close to the northern auroral region, the long- 
wavelength emission marks it out as unique when compared with auroral features. It is 
also unique compared to the transient bright flashes described in Giles et al. (2020), which 
were all dominated by Hə emission, like the aurora. Despite these different spectral char- 
acteristics, the bright spot does clearly originate from within Jupiter’s atmosphere; Fig- 
ure 1 shows that it is located far from the limb of the planet so there is no possibility 
of confusion with a star close to Jupiter, and the bright spot spectrum shows CH4 ab- 
sorption from Jupiter’s atmosphere (see Section 3.2) so it cannot originate from an ob- 


ject located between the spacecraft and the planet. 


We conducted a search for other such events, by filtering for occasions when there 
was a sharp increase and then decrease in the number of long-wavelength photon counts 
recorded, but no similar events were found. The spot is analyzed further in Section 3 and 


the results are discussed in Section 4. 


3 Results 
3.1 Spatial extent and duration 


Because of the way in which the UVS maps are built up as the spacecraft spins, 
Figure 1 contains information about both the spatial extent and the duration of the bright 
spot. By fitting a two-dimensional Gaussian to the total photon counts, we find that the 


spot has a FWHM of 0.29° in the along slit direction and 0.23° in the across slit direc- 


tion. This is consistent with the shape that we would expect for a point source that has 

a constant brightness on the 17-ms timescale it takes the slit to pass over the source (Greathouse 
et al., 2013). Based on the distance of the spacecraft at the time of the observations (80,000 
km), the upper limit for the spot diameter is 400 km. 


Because the shape of the bright spot is consistent with the FWHM of a point source, 
the bright spot must also be approximately stationary in latitude and longitude over the 
17-ms duration of the observation. In order to observe motion, the spot would have had 
to move >0.1° in 17 ms, equating to a horizontal velocity of >8000 kms"!. Since we do 


not observe any elongation, any horizontal movement must be below this speed. 


While the brightness appears approximately constant on the timescale it takes the 
slit to pass over the source, the bright spot was only observed during a single spin, so 
it must be transient on slightly longer timescales. UVS data was recorded from the same 
latitude and longitude two spins (60 s) earlier and three spins (90 s) later, and the bright 
spot was not present in either of these maps. This places an upper bound of 150 s on 


the duration of the bright spot, giving a duration of between 17 ms and 150 s. 


Because the bright spot was observed on the night side of the planet, it is not pos- 
sible to use UVS data to search for any local atmospheric changes in its aftermath, such 
as changes in the aerosol abundance. Any such changes could therefore have lasted longer 


than 150 s. 


3.2 Spectral analysis 


The bright spot spectrum was calculated by summing all photon detections within 
a 1°x1° box centered on the bright spot. The spectrum is shown by the black data points 
in Figure 2 and is presented in terms of spectral irradiance measured at the spacecraft. 
As suggested earlier by the color of the bright spot in Figure 1, the bright spot emission 
has an unusual spectral shape compared to both the auroral emission and the TLE emis- 
sion described in Giles et al. (2020). The auroral and TLE spectra are dominated by Hə 
emission, and the Hə Lyman band can be recognized by a double peak at 160 nm. Be- 
yond 170 nm, the Hz emission drops to low levels. In contrast, the spectrum shown in 
black in Figure 2 increases smoothly between 150 and 190 nm, suggesting blackbody emis- 
sion. The rapid drop-off at wavelengths shorter than 140 nm suggests strong atmospheric 


absorption by CH4. 


In order to model a blackbody spectrum attenuated by atmospheric absorption, 
we used the atmospheric composition Model C from Moses et al. (2005). In this model 


atmosphere, the CH, homopause is located at a pressure of 2.9x10~* mbar, or 360 km 
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Figure 2. The bright spot spectrum (black) compared with various model spectra. (a) The 


colored lines show the shape of blackbodies with different temperatures. There is no atmospheric 
absorption. (b) The colored lines show the spectral shape of a 9000 K blackbody placed at dif- 
ferent altitudes and therefore with different amounts of atmospheric absorption. (c) The colored 


lines show two different best-fit model spectra. 


above the 1 bar level. For a given altitude and emission angle, we used the atmospheric 
composition model to calculate the column density for each atmospheric gas. The three 
stratospheric gases that contribute significantly to absorption in the 125-210 nm range 
are CH4, CoH, and C2H¢g. Absorption cross-sections for these three gases were obtained 
from Chen and Wu (2004), Lee et al. (2001), Smith et al. (1991) and Bénilan et al. (2000). 
The atmospheric transmission was then calculated using the column density for each gas 
and the absorption cross-sections, and this was multiplied by a blackbody spectrum to 
obtain a modeled top-of-atmosphere spectrum. The colored lines in Figures 2(a) and 2(b) 
show these modeled spectra for a range of blackbody temperatures and emission alti- 


tudes. 


Figure 2(a) shows modeled spectra for blackbodies with five different temperatures. 
For these spectra, no atmospheric absorption is included and the spectra have been scaled 
to match the average irradiance of the bright spot spectrum. Within the 125-210 nm 
spectral range, increasing the temperature from 7000 K to 12,000 K flattens the spec- 
trum. At shorter wavelengths, the data diverges sharply from all of the blackbody spec- 
tra due to the lack of atmospheric absorption. We note that the data also diverges from 
the blackbody spectra at ~195 nm. As demonstrated by the larger error bars, the UVS 
sensitivity is much lower at these longer wavelengths than at shorter wavelengths; by 195 
nm, the instrument’s effective area is an order of magnitude smaller than at 160 nm (Hue 
et al., 2018). We therefore attribute this apparent small dip in the irradiance to limi- 


tations in the radiometric calibration. 


Figure 2(b) shows modeled spectra obtained using different emission altitudes (i.e. 
different amounts of atmospheric absorption). Each spectrum uses the same 9000 K black- 
body as the source emission and the altitude is defined as the height above the 1-bar level. 
The red lines in Figures 2(a) and 2(b) are identical, as there is negligible absorption above 
an altitude of 400 km. As the source is moved deeper in the atmosphere, the amount of 
absorption increases. By 250 km there is significant CH, absorption shortwards of 140 
nm. By 200 km, there is also a significant C2H2 feature at 152 nm. By 150 km, the at- 
mosphere is opaque shortwards of 155 nm due to the combination of C2H2, C2H¢ and 


CH4. 


When the blackbody temperature, emission altitude and a scaling factor were all 


allowed to vary together, the best fit was obtained with a temperature of 9600+600 K 


and an altitude of 225+5 km above the 1-bar level (assuming an error of 20% on the abun- 
dances of the absorbers in the model). This best-fit spectrum is shown by the red line 
in Figure 2(c). While this modeled spectrum provides a good fit at most wavelengths, 


the C2Hə absorption feature at 152 nm is clearly too strong; an altitude of 225 km pro- 


vides a good fit for the CH, and C2H¢ column densities, but the C2H2 column density 
(2.1x10'© cm) is too high. Instead, if we allow the C2H»2 column density to vary in- 
dependently of the other gases, we retrieve a value of 1.1 x 1016 cm, approximately 
half of the value at 225 km. The spectrum obtained with this value is shown by the blue 
line in Figure 2(c). Compared to the Moses et al. (2005) model atmosphere, this sug- 
gests that the atmosphere near the bright spot has lower C2H2 volume mixing ratio. This 
is unsurprising as the Moses et al. (2005) model describes equatorial latitudes, and Nixon 
et al. (2007) used Cassini CIRS observations to show that the C2H2 volume mixing ra- 
tio decreases at higher latitudes. The blue line in Figure 2(c) is therefore consistent with 


an altitude of 225 km at a latitude of 53°N. 


4 Discussion 


The bright spot described in Section 3 has the following properties: 


(i) Its duration is longer than 17 ms and less than 150 s. 
(ii) The shape is consistent with a point source, which equates to a maximum diam- 
eter of 400 km. 
(iii) Its spectral shape is consistent with a blackbody of temperature 9600 K located 
at an altitude of 225 km above the 1-bar level (a pressure of 0.04 mbar). 


In Giles et al. (2020), we presented several instances of transient bright flashes ob- 
served in the Juno UVS data and we concluded that these were consistent with TLEs 
known as elves, sprites or sprite halos, upper atmospheric responses to tropospheric light- 
ning. The bright spot that we present in this paper differs from the previous bright flashes 
in two main ways, highlighted in Figure 3. Firstly, it has a significantly longer duration. 
The TLEs had a duration of ~1.4 ms, which is short enough that UVS could observe their 
decay during the time it took the instrument slit to pass across the point source. In con- 
trast, this bright spot is approximately constant as the slit passes across it, giving it a 
minimum duration of ~17 ms. This difference in duration is shown in Figures 3(a) and 
(b); the Gaussian shape in (a) is governed by the width of the slit, while the rapid ex- 
ponential decay of (b) takes place on a much shorter timescale. Secondly, the spectra 
of the TLEs were dominated by Hz Lyman band emission, giving them a very similar 
spectral shape to auroral emission. In contrast, this bright spot has a blackbody-like emis- 
sion spectrum. The spectra are compared in Figures 3(c) and (d). As discussed in Sec- 
tion 3, the shape of (c) is consistent with blackbody emission, while the double peak at 


160 nm seen in (d) is indicative of Hz emission. 
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Figure 3. <A comparison of the time series and spectra of the bolide (this paper) and the 

TLEs (Giles et al., 2020). (a) and (b) show the number of observed UV photons as a function 
of time, relative to the peak event time. The time series shown in (b) is for the TLE observed 
on April 10 2020 at 17:24:35 UTC. (c) shows the same spectrum presented in Figure 2 and (d) 


shows the averaged TLE spectrum from Giles et al. (2020). 
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Because of both the duration and the spectral shape, the bright spot described in 
this paper is unlikely to be a TLE. Lightning and associated TLEs do not have black- 
body spectra; on earth, they have emission spectra that are dominated by nitrogen and 
oxygen lines (Walker & Christian, 2017; Rodger, 1999) and on Jupiter their emission spec- 
tra are dominated by hydrogen (Borucki et al., 1996; Yair et al., 2009). Elves have sub- 
millisecond timescales, and while sprites can last for several tens of milliseconds, the brighter 
ones are often shorter (Lyons, 2006). There are other longer-lasting TLEs, such as blue 
jets, but these would still be Hz dominated and they emerge directly from the thunder- 
storm anvils (Rodger, 1999), so we would expect them to be deeper in the atmosphere. 
The shape of the spectrum also rules out the possibility that this is a transient auroral 


event, as auroral emission is Hy dominated. 


Given the spectrum and the duration, we suggest that the bright spot we observe 
may be a bolide/fireball in Jupiter’s atmosphere. Borovička and Spurný (1996) studied 
the lightcurves and spectra of two of the brightest bolides observed in the Earth’s at- 
mosphere and found that the spectra were a mixture of blackbody continuum and emis- 
sion lines. The emission lines come from ablated meteoroid material and the blackbody 
continuum comes from thermal bremsstrahlung (free-free) emission. In the beginning phase 
of the bolide, Borovička and Spurný (1996) found that the line emission dominated, but 
during the brightest part of the lightcurve, the blackbody continuum dominated. A tem- 
perature of 9600 K is approximately consistent with the Planck temperatures measured 
in Earth bolides. Borovicka (1994) found that fireball spectra in the Earth’s atmosphere 
exhibit two distinct characteristic temperatures: a primary component with a brightness 
temperature of 4000 K and a secondary component with a brightness temperature of 10,000 
K. The lower temperature is thought to be from the thermal radiation of the meteoroid 
itself, while the second part is associated with the shock wave. The greater the speed 
of the meteoroid, the more dominant this second part is. We do not observe any move- 
ment of the bright spot during our 17-ms observation timescale, and this is consistent 
with a typical impact velocity of tens of kms! (Crawford et al., 1994); a horizontal ve- 


locity of >8000 kms"! would be required to be detectable by UVS. 


Our observations are somewhat similar to the observations made by the UVS in- 
strument on the Galileo spacecraft when comet Shoemaker-Levy 9 collided with Jupiter 
in 1994 (Hord et al., 1995). They observed a brightening during one spatial scan only 
(no brightening in scans obtained 5 s beforehand or 5 s afterwards) and when combined 
with simultaneous observations from the Photopolarimeter Radiometer, they concluded 
that the brightness temperature was 78007 200 K. Our observations are also similar to 


several observations made by amateur observers of bolides on Jupiter (Hueso et al., 2018). 


[i= 


These observations were recorded in the visible and the observed bright flashes lasted 
1-2 s in each case. One of these events was simultaneously recorded in two different fil- 
ters, and fitting a blackbody to these two measurements gave a temperature of 6500- 


8500 K (Hueso et al., 2013). 


Hueso et al. (2010) used the Earth-based observation of a Jovian bolide to estimate 
the size of the impacting object and we follow their analysis approach here. If we cor- 
rect for atmospheric absorption, extend the blackbody shape to all wavelengths and in- 
tegrate, the total irradiance observed at the spacecraft would be 1.2—2.3x107!° Wem”? 
(taking into account the 600 K error on the brightness temperature). Assuming isotropic 
radiation, the total power emitted is 1.0—1.9x10!! W. As we only observe the bright 
flash for a very short period of time as the spacecraft spins, we do not know the total 
duration of the emission and this adds significant uncertainty to the calculation. Here, 
we assume a duration of 1-2 s, based on Hueso et al. (2018). There is additional uncer- 
tainty from the fact that we do not know which stage of the light curve we are observ- 
ing - it could be the peak of the emission, or the peak could occur shortly before/after 
our observations. Including a 50% uncertainty for this factor, we obtain a total optical 


energy of 5 x 10!°-6 x 101! J. 


Brown et al. (2002) found that the efficiency, u, with which the kinetic energy of 


an impactor is converted into optical energy can be empirically described by 


p =0.121 x E915 (1) 


where Ep is the optical energy in terms of kilotons of TNT (1 kton = 4.185x101? 
J). For our observed total optical energy range, the efficiency is 7-10% and the kinetic 
energy of the impactor is therefore 7 x 1011- 6 x 10!” J. Assuming an impactor veloc- 
ity of 60 kms", the velocity at which fragments of Shoemaker-Levy 9 collided with Jupiter (Crawford 
et al., 1994), and a velocity error of 20% leads to a mass estimate of 250-5000 kg. For 


densities of of 250-2000 kgm’, this equates to a diameter of 1-4 m. 


By considering the total power emitted and the blackbody temperature, we can use 
the Stefan-Boltzmann law to calculate the area of the emitting region. Assuming an emis- 
sivity of 1, the effective diameter of the emitting region is 9 m. This is consistent with 
an impactor diameter of 1-4 m; Borovička and Spurný (1996) found that the maximum 
diameter of the radiating region is on the order of 10 times larger than the initial diam- 


eter of the body. 


== 


A mass estimate of 250-5000 kg seems broadly consistent with the altitude at which 
we observe the bright flash: 225 km above the l-bar level, or a pressure of 0.04 mbar. 
Numerical impact simulations have been used to estimate the altitudes of bolides in Jupiter’s 
atmosphere, and have found that impactors with masses of ~ 1 x 10° kg reach their 
peak brightness at 60-150 km above the 1 bar level (Hueso et al., 2013; Sankar et al., 
2020). Smaller impactors, such as our observation, would not penetrate as deeply into 
the atmosphere. Borovička and Spurný (1996) studied the entry of a 5000-kg impactor 
in the Earth’s atmosphere, equivalent to the upper end of our mass estimate. They found 
that the peak brightness occurred at an altitude of 67 km, which corresponds to a pres- 
sure of 0.07 mbar (ISO, 1975). This is slightly deeper in the atmosphere than our ob- 
served bright flash pressure level of 0.04 mbar, which is consistent with 5000 kg being 


the upper bound of our mass estimate. 


By considering all observations made with Juno UVS over the first 27 perijoves of 
the mission, we find that UVS obtained a total effective coverage of 8.2 x 10'8 km?s. 
In that time, there was a single bolide observation. There are clear limitations in cal- 
culating an impact flux rate from a single observation, but a maximum likelihood esti- 
mation calculation leads to an impact rate of ~24,000 per year (See Text S1 in the Sup- 
porting Information). By considering how impact rates vary as a function of meteoroid 
size, we can compare this impact rate with the rates deduced from previous studies. We 
find that our rate is lower than the flux rate estimated by Cook and Duxbury (1981) and 
higher than the flux rate of Hueso et al. (2018) (See Text S2 in the Supporting Infor- 
mation). This comparison is dependent on the assumed mass of our observed meteor, 
and the lower limit of 250 kg is significantly more consistent with Hueso et al. (2018) than 
the upper limit of 5000 kg. We note that there are considerable uncertainties in our es- 


timated rate, as we have thus far only observed a single bolide with Juno UVS. 


5 Conclusions 


Juno UVS observations recorded transient blackbody emission from a point source 
in Jupiter’s atmosphere. Spectral modelling showed that the emission is consistent with 
a 9600 K source located 225 km above the 1-bar level, and the emission lasted between 
17 ms and 150 s. The blackbody nature of the spectrum, the temperature and the du- 
ration of the emission are all consistent with a bolide in Jupiter’s atmosphere. Based on 
the energy emitted, we estimate that the impactor had a mass of 250-5000 kg, which would 
correspond to a diameter of 1-4 m. This impactor size is larger than the small fireball 
observed by Cook and Duxbury (1981) and smaller than the superbolides described by 


(Hueso et al., 2018). We estimate an impact flux rate of 24,000 per year, for masses of 


=[3= 


250-5000 kg or greater; if we use the lower limit of our mass estimate, this is consistent 
with the rate determined by Hueso et al. (2018). However, we note that there are large 
uncertainties associated with this value due to the fact that we have thus far only ob- 
served one bolide with Juno UVS. As the Juno mission continues, our integration time 
on the planet will continue to steadily increase and our estimate of the impact flux rate 


will improve. 


Acknowledgements 


We are grateful to NASA and contributing institutions, which have made the Juno 
mission possible. This work was funded by NASA’s New Frontiers Program for Juno via 
contract with the Southwest Research Institute. B.B. is a Research Associate of the Fonds 


de la Recherche Scientifique - FNRS. 


Data Availability Statement 


The Juno UVS data used in this paper are archived in NASA’s Planetary Data Sys- 
tem Atmospheres Node: https://pds-atmospheres.nmsu.edu/PDS/data/jnouvs_3001 (Trantham, 
2014). The data used to produce the figures in this paper are available in Giles (2021). 


References 


Adriani, A., Coradini, A., Filacchione, G., Lunine, J. I., Bini, A., Pasqui, C., ... 
others (2008). JIRAM, the image spectrometer in the near infrared on board 
the Juno mission to Jupiter. Astrobiology, 8(3), 613-622. 

Asphaug, E., & Benz, W. (1996). Size, density, and structure of comet Shoemaker- 
Levy 9 inferred from the physics of tidal breakup. Icarus, 121(2), 225-248. 
Bénilan, Y., Smith, N., Jolly, A., & Raulin, F. (2000). The long wavelength range 
temperature variations of the mid-UV acetylene absorption coefficient. Plane- 

tary and Space Science, 48 (5), 463-471. 

Benmahi, B., Cavalié, T., Dobrijevic, M., Biver, N., Bermudez-Diaz, K., Sandqvist, 
A., ... others (2020). Monitoring of the evolution of H20 vapor in the strato- 
sphere of Jupiter over an 18-yr period with the Odin space telescope. Astron- 
omy & Astrophysics, 641, A140. 

Bland, P. A., & Artemieva, N. A. (2006). The rate of small impacts on Earth. Mete- 
oritics & Planetary Science, 41 (4), 607-631. 

Bolton, S. J., Lunine, J., Stevenson, D., Connerney, J. E. P., Levin, S., Owen, T. C., 

. others (2017). The Juno Mission. Space Science Reviews, 213(1-4), 5-37. 

Bonfond, B., Gladstone, G. R., Grodent, D., Greathouse, T. K., Versteeg, M. H., 


[A= 


Hue, V.,... others (2017). Morphology of the UV aurorae Jupiter during 
Juno’s first perijove observations. Geophysical Research Letters, 44 (10), 4463- 
4471. 

Borovička, J. (1994). Two components in meteor spectra. Planetary and Space Sci- 
ence, 42 (2), 145-150. 

Borovička, J., & Spurný, P. (1996). Radiation study of two very bright terrestrial 
bolides and an application to the comet S-L 9 collision with Jupiter. Icarus, 
121 (2), 484-510. 

Borucki, W. J., McKay, C. P., Jebens, D., Lakkaraju, H. S., & Vanajakshi, C. T. 
(1996). Spectral irradiance measurements of simulated lightning in planetary 
atmospheres. Icarus, 123 (2), 336-344. 

Brown, P., Spalding, R. E., ReVelle, D. O., Tagliaferri, E., & Worden, S. P. (2002). 
The flux of small near-Earth objects colliding with the Earth. Nature, 
420 (6913), 294-296. 

Cavalié, T., Feuchtgruber, H., Lellouch, E., de Val-Borro, M., Jarchow, C., Moreno, 
R., ... others (2013). Spatial distribution of water in the stratosphere of 
Jupiter from Herschel HIFI and PACS observations. Astronomy & Astro- 
physics, 553, A21. 

Chen, F., & Wu, C. R. (2004). Temperature-dependent photoabsorption cross sec- 
tions in the VUV-UV region. I. Methane and ethane. Journal of Quantitative 
Spectroscopy and Radiative Transfer, 85(2), 195-209. 

Cook, A. F., & Duxbury, T. C. (1981). A fireball in Jupiter’s atmosphere. Journal of 
Geophysical Research: Space Physics, 86(A10), 8815-8817. 

Crawford, D. A., Boslough, M. B., Trucano, T. G., & Robinson, A. C. (1994). The 
impact of comet Shoemaker-Levy 9 on Jupiter. Shock Waves, 4(1), 47-50. 

Eliason, S. R. (1993). Maximum likelihood estimation: Logic and practice. Sage. 

Giles, R. S. (2021). Detection of a bolide in Jupiter’s atmosphere with Juno UVS. 
Mendeley Data, V1. doi: 10.17632/nhc5r85xv3.1 

Giles, R. S., Greathouse, T. K., Bonfond, B., Gladstone, G. R., Kammer, J. A., Hue, 
V.,... Levin, S.M. (2020). Possible Transient Luminous Events observed 
in Jupiter’s upper atmosphere. Journal of Geophysical Research: Planets, 
125(11), e2020JE006659. 

Gladstone, G. R., Persyn, S. C., Eterno, J. S., Walther, B. C., Slater, D. C., Davis, 
M. W., ... others (2017). The ultraviolet spectrograph on NASA’s Juno 
mission. Space Science Reviews, 213(1-4), 447-473. 

Gladstone, G. R., Versteeg, M. H., Greathouse, T. K., Hue, V., Davis, M. W., 
Gérard, J.-C., ... others (2017). Juno-UVS approach observations of Jupiter’s 


jhe 


auroras. Geophysical Research Letters, 44(15), 7668-7675. 

Greathouse, T. K., Gladstone, G. R., Davis, M. W., Slater, D. C., Versteeg, M. H., 
Persson, K. B., ... Eterno, J. S. (2013). Performance results from in-flight 
commissioning of the Juno Ultraviolet Spectrograph (Juno-UVS). In UV, X- 
Ray, and Gamma-Ray Space Instrumentation for Astronomy XVIII (Vol. 8859, 
p. 88590T). 

Hansen, C. J., Caplinger, M. A., Ingersoll, A., Ravine, M. A., Jensen, E., Bolton, 

S., & Orton, G. (2017). Junocam: Juno’s outreach camera. Space Science 
Reviews, 213 (1-4), 475-506. 

Hord, C. W., Pryor, W. R., Stewart, A. I. F., Simmons, K. E., Gebben, J. J., Barth, 
C. A.,... others (1995). Direct observations of the comet Shoemaker-Levy 
9 fragment G impact by Galileo UVS. Geophysical Research Letters, 22 (12), 
1565-1568. 

Hue, V., Kammer, J., Gladstone, G. R., Greathouse, T. K., Davis, M. W., Bonfond, 
B.,... others (2018). In-flight characterization and calibration of the Juno- 
Ultraviolet Spectrograph (Juno-UVS). In Space telescopes and instrumentation 
2018: Ultraviolet to gamma ray (Vol. 10699, p. 1069931). 

Hueso, R., Delcroix, M., Sánchez-Lavega, A., Pedranghelu, S., Kernbauer, G., McK- 
eon, J., ... others (2018). Small impacts on the giant planet Jupiter. Astron- 
omy & Astrophysics, 617, A68. 

Hueso, R., Pérez-Hoyos, S., Sánchez-Lavega, A., Wesley, A., Hall, G., Go, C., ... 
others (2013). Impact flux on Jupiter: From superbolides to large-scale colli- 
sions. Astronomy & Astrophysics, 560, A55. 

Hueso, R., Wesley, A., Go, C., Pérez-Hoyos, S., Wong, M., Fletcher, L., ... oth- 
ers (2010). First Earth-based detection of a superbolide on Jupiter. The 
Astrophysical Journal Letters, 721(2), L129. 

ISO. (1975). ISO 2533:1975 - Standard Atmosphere. International Organization for 
Standardization. 

Lee, A. Y., Yung, Y. L., Cheng, B.-M., Bahou, M., Chung, C.-Y., & Lee, Y.-P. 
(2001). Enhancement of deuterated ethane on Jupiter. The Astrophysical 
Journal Letters, 551(1), L93. 

Lellouch, E., Bézard, B., Strobel, D. F., Bjoraker, G. L., Flasar, F. M., & Romani, 
P. N. (2006). On the HCN and CO2 abundance and distribution in Jupiter’s 
stratosphere. Icarus, 184(2), 478-497. 

Lyons, W. A. (2006). The meteorology of Transient Luminous Events — an introduc- 
tion and overview. In Sprites, Elves and Intense Lightning Discharges (pp. 19- 
56). Springer. 


—16- 


Moses, J. I., Fouchet, T., Bézard, B., Gladstone, G. R., Lellouch, E., & Feucht- 
gruber, H. (2005). Photochemistry and diffusion in Jupiter’s stratosphere: 
constraints from ISO observations and comparisons with other giant planets. 
Journal of Geophysical Research: Planets, 110(E8), E08001. 

Nixon, C. A., Achterberg, R. K., Conrath, B. J., Irwin, P. G. J., Teanby, N. A., 
Fouchet, T., ... others (2007). Meridional variations of C2H2 and C2Hg in 
Jupiter’s atmosphere from Cassini CIRS infrared spectra. Icarus, 188 (1), 
47-71. 

Poppe, A. R. (2016). An improved model for interplanetary dust fluxes in the outer 
Solar System. Icarus, 264, 369-386. 

Rodger, C. J. (1999). Red sprites, upward lightning, and VLF perturbations. Re- 
views of Geophysics, 87(3), 317-336. 

Sdnchez-Lavega, A., Gómez, J. M., Rojas, J. F., Acarreta, J. R., Lecacheux, J., Co- 
las, F., ... Arregui, J. (1998). Long-term evolution of comet SL-9 impact 
features: July 1994-September 1996. Icarus, 131(2), 341-357. 

Sánchez-Lavega, A., Orton, G. S., Hueso, R., Pérez-Hoyos, S., Fletcher, L. N., 


García-Melendo, E.,... others (2011). Long-term evolution of the aerosol 
debris cloud produced by the 2009 impact on Jupiter. Icarus, 214 (2), 462- 
476. 


Sánchez-Lavega, A., Wesley, A., Orton, G., Hueso, R., Perez-Hoyos, S., Fletcher, L., 
. others (2010). The impact of a large object on Jupiter in 2009 July. The 

Astrophysical Journal Letters, 715(2), L155. 

Sankar, R., Palotai, C., Hueso, R., Delcroix, M., Chappel, E., & Sanchez-Lavega, 
A. (2020). Fragmentation modelling of the 2019 August impact on Jupiter. 
Monthly Notices of the Royal Astronomical Society, 498 (4), 4622-4630. 

Smith, P. L., Yoshino, K., Parkinson, W., Ito, K., & Stark, G. (1991). High- 
resolution, VUV (147-201 nm) photoabsorption cross sections for C2Hə at 195 
and 295 K. Journal of Geophysical Research: Planets, 96 (E2), 17529-17533. 

Sremčević, M., Krivov, A. V., Krüger, H., & Spahn, F. (2005). Impact-generated 
dust clouds around planetary satellites: model versus Galileo data. Planetary 
and Space Science, 58 (6), 625-641. 

Trantham, B. (2014). JNO-J-UVS-3-RDR-V1.0. NASA Planetary Data System. 

Walker, T. D., & Christian, H. J. (2017). Triggered lightning spectroscopy: Part 
1. A qualitative analysis. Journal of Geophysical Research: Atmospheres, 
122 (15), 8000-8011. 

Yair, Y., Takahashi, Y., Yaniv, R., Ebert, U., & Goto, Y. (2009). A study of the 
possibility of sprites in the atmospheres of other planets. Journal of Geophysi- 


=[7= 


cal Research: Planets, 114(E9), E09002. 


-jg 


Supporting Information for “Detection of a bolide in 
Jupiter’s atmosphere with Juno UVS” 


Contents of this file 


1. Text S1: Impact rate 


2. Text S2: Comparison of impact rate with previous studies 


Introduction 


This supporting information contains details about the impact flux rate calcula- 


tion and the comparison to previous impact flux rate estimates. 
Text S1: Impact rate 


By multiplying the time spent observing Jupiter by the footprint of the UVS slit 
on the planet, we find that over 27 perijoves, UVS obtained coverage of 9.3x101! km?s. 
However, a given point on the planet is only observed for 17 ms as the instrument slit 
passes over it, and the bright flash caused by a meteor typically lasts 1-2 s. Therefore, 
the effective coverage considered in the search for bright flashes must be 1.5/0.017 = 88 
times larger, giving a total effective coverage of 8.2x10!° km?s. Using the surface area 
of the planet, this is the equivalent of observing the entirety of Jupiter for 1300 seconds 


(although we note that our observations are not evenly distributed across the planet). 


During that time period, we observed one bolide. It is not possible to calculate an 
accurate estimate of the average occurrence rate from a single event. However, we can 
use maximum likelihood estimation (Eliason, 1993) to find the occurrence rate that is 
most likely to produce a single bolide observation in our observation time period. If bolides 
occur independently of each other, their occurrence can be described by a Poisson dis- 
tribution. In this case, the probability of a single event occurring in time T is given by 
(rT)e—"?, where r is the occurrence rate. This probability is a maximum when r = 1/T. 


In this case, T is 1300 seconds, so r is 1 per 1300 seconds, or ~24,000 per year. 


Our observation was made from a distance of 80,000 km and we estimate that the 
bolide we observed was due to an impactor of mass 250-5000 kg. While we would cer- 
tainly be able to detect smaller impactors when the spacecraft is even closer to the planet, 
much of the observing time is when the spacecraft is further away. It is therefore diffi- 
cult to associate the impact flux with a specific mass threshold; in the following discus- 


sion, we use both the 250 kg and 5000 kg values. 
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Text S2: Comparison of impact rate with previous studies 


Bland and Artemieva (2006) compiled a range of observational studies in order to 
produce a graph showing the rate of impact of meteoroids in the Earth’s upper atmo- 
sphere as a function of meteoroid size. There is a log-log relationship between impact 


rate and size. For meteoroids in the mass range 3-10!° kg, the relationship is 
log N = —0.926 log m + 4.739 (2) 


where N is the number of impacts per year of mass m or greater. Using this equation, 

the Earth undergoes ~330 impacts of size >250 kg and ~20 impacts of size >5000 kg 
each year. However, given the difference in orbit, mass and area, we clearly expect flux 
rates to be different on Jupiter. There have been two previous estimates of small impactor 
flux rates on Jupiter, one based on even smaller impactors (Cook & Duxbury, 1981), and 
one based on slightly larger impactors (Hueso et al., 2018). If we assume that the dis- 
tribution of meteoroid sizes is the same for Jupiter as it is for the Earth, we can use Equa- 


tion 2 to compare our flux rate with these previous studies. 


Cook and Duxbury (1981) used the camera on Voyager 1 to search for fireballs and 
detected one such event. Based on this, they estimated an impact rate of 1.6 x 10719 
stkm? (3 x 108 impacts per year) for impactors of mass >2.8 kg. Using Equation 2, 
we would expect >250 kg impactors to be 64 times less frequent than >2.8 kg impactors, 
i.e. the results of Cook and Duxbury (1981) suggest a rate of 5x10° impacts per year 
for impactors >250 kg. Similarly, we would expect >5000 kg impactors to be 1000 times 
less frequent than >2.8 kg impactors, i.e. the results of Cook and Duxbury (1981) sug- 
gest a rate of 3x10° impacts per year for impactors >5000 kg. Our estimated flux rate 


of 24,000 impacts per year is approximately 10-200 times lower than these values. 


Hueso et al. (2018) combine the results of several amateur observations over mul- 
tiple years to estimate a impact rate for Jupiter for impactors in the 1 x 10°—1 x 10® 
kg range. Amateur observers have observed five fireballs in Jupiter’s atmosphere in 2010- 
2017, and by estimating the effective time during which the amateur observing commu- 
nity takes measurements, Hueso et al. (2018) estimate an impact rate of 10-65 events 
per year. Following the mass distribution given in Equation 2, this corresponds to a rate 
of 2,600-17,000 impacts per year for masses >250 kg and a rate of 160—-1,000 impacts 
per year for masses >5000 kg. Our estimated flux rate of 24,000 impacts per year is higher 
than both of these values, but our lower mass estimate of 250 kg produces an impact rate 


that is significantly more consistent with Hueso et al. (2018). 


The previous two paragraphs compare three impact rate estimates: 3 x 108 per 


year for impactors >2.8 kg (Cook & Duxbury, 1981), 24,000 per year for impactors >250- 


~20- 


5000 kg (this paper) and 10-65 per year for impactors >1x10° kg (Hueso et al., 2018). 
If the mass distribution of impactors follows Equation 2, then our measured rate is lower 
than the value expected from Cook and Duxbury (1981) and higher than the value ex- 
pected from Hueso et al. (2018). However, we note that all three rates have a significant 
amount of uncertainty as they are based on a very small number of impact observations: 
one each in the case of this paper and Cook and Duxbury (1981) and five in the case of 
Hueso et al. (2018). While Hueso et al. (2018) has a larger number of impacts, the ef- 


fective duration of the observations is uncertain. 
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